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Abstract. - We perform an orthogonal basis tight binding fit to an LAPW calculation of 
paramagnetic Na^CoCh for several dopings. The optimal position of the apical oxygen at 
each doping is resolved, revealing a non-trivial dependence of the band structure and Fermi 
surface on oxygen height. We find that the small e g i hole pockets are preserved throughout all 
investigated dopings and discuss some possible reasons for the lack of experimental evidence 
for these Fermi sheets. 



The compound Na^CoC^ has been synthesized for a wide range of Na contents, 0.3 < x 
< 0.9, and exhibits a variety of unusual phenomena that are strongly dependent on doping 
level. Study of Nao.5Co02 for its surprisingly good thermoelectric properties began sev- 
eral years ago [1] and more recently, attention has been focused on the hydrated compound 
Nao.35CoO2T.3H2O which undergoes a superconducting phase transition at T ~ 4.5K [2,3]. 

The behavior and characteristics of Naa:Co02 do not vary smoothly along the spectrum 
of Na contents, but instead split very generally into two separate metallic phases separated 
at x=0.5 by an insulating phase [4]. The low Na content region at x <0.5 is characterized by 
conventional metallic behavior [5] and a nearly temperature independent susceptibility [6, 8] 
indicating a Pauli paramagnetic precursor to the superconducting state. The linear specific 
heat coefficient 7 is measured to be ~ 12-16 m T " ; ^y which shows only light mass renormalization 
[6,7]. In contrast, the high Na region at x > 0.5, though also metallic, exhibits a Curie- Weiss- 
like susceptibility [9,10,12]. A higher (compared to the previous metallic phase) measured 
7 of 25-30 m ™^. a indicates that mass renormalization in this region is substantial [9, 12, 15]. 
Interestingly, the Kadawaki- Woods ratio, a comparison of 7 to the quadratic coefficient of 
resisitivity, which is often used as a measure of electron-electron interaction, is reported to 
be the largest ever measured and considerably field dependent [16]. This may evidence a 
nearby quantum critical point and indeed, fisR and transport measurements show a magnetic 
transition in single crystals at T=22K for dopings of x= 0.7 and x= 0.75, usually interpreted 
as a spin density wave [10, 13, 14]. Strong ferromagnetic spin fluctuations have been predicted 
based on LDA calculations [11] and have been detected by neutron scattering [17] when x= 
0.7, but with slightly higher electron count at x=0.85, a bulk anti-ferromagnetic transition has 
been reported [18] , suggesting the presence of more than one type of magnetic interaction. The 
insulating phase is very narrowly centered at x=0.5 and recent reports show both charge and 
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magnetic order [4], though earlier investigations, most notably work by Terasaki et al [1,19,20], 
found this system to be metallic. Measurements of 7 ~ 40-56 m ";^a are higher than for 
any other doping [19, 20] and NMR measurements indicate Curie- Weiss behavior and charge 
separation [21]. 

Here we present a highly accurate tight-binding parameterization of Na K Co02, obtained 
for a range of apical oxygen heights at each of three dopings: x=0.3, 0.5, and 0.7, meant to be 
representative of the three general phases. We anticipate that this Hamiltonian will be valuable 
as accurate, material-specific input to model systems that can incorporate the important 
effects of correlation. We point out that, since the system is multi-band in nature, a strictly 
two-dimensional triangular lattice model is insufficient to fully account for all dispersions. In 
contrast to earlier tight-binding models [22,23], we include both interlayer couplings and all 
five Co-d orbitals as a basis. The nearly perfect reproduction of LDA Fermi level crossings 
within our tight-binding scheme allows a detailed and accurate examination of doping and 
structural Fermi surface dependencies. We show that the band structure is strongly coupled 
to the apical oxygen height which itself depends on the amount of charge donated to the C0O2 
plane by Na ions. A rigid band picture is insufficient to describe the evolution from small 
x where the band structure is nearly two dimensional to large x where non-negligible inter- 
planar coupling causes heavy distortion of one Fermi sheet. An important effect of increased 
three-dimensionality in the high Na compounds is the preservation of the small hole pockets 
at the K points of the Brillouin zone (BZ). We discuss possible reasons that these pockets are 
absent from ARPES measurements [28,29]. 

Through a series of calculations with an LAPW implementation (Wien2k) [24] of DFT 
formalism, we varied the Na level, x, of Na^OoC^, holding the total volume constant at 
509.553 a.u. 3 . Though important effects may derive from the specific location of Na ions, 
particularly in the insulating compound at x—0.5, where Na and charge order are thought to 
cooperate, we treat Na in the virtual crystal approximation (VCA) and concentrate only on 
effects stemming directly from variation of electron number. Specifically, we treat every Na 
site as fully occupied, but with a fictitious ion of atomic charge 10.x, rather than a true Na 
atom of charge 11. This approximation is most valid when the Na states are far from the 
Fermi energy, as they are Na^CoC^ for all x, and in fact, comparison of VCA and Na ordered 
calculations show little difference [25]. The apical oxygen height was relaxed for a number of 
dopings and shown to increase non-linearly with Na addition. The optimal O position as a 
fraction of the lattice constant (a= 5.366824 a.u.) for x = 0.3, 0.5 and 0.7 was found to be 
0.0811, 0.0834, and 0.0864 respectively. The six Co-derived bands near the Fermi energy in 
each of these first-principles calculations form the basis for our tight-binding fit. 

We fit the LAPW eigenvalues using the NRL tight- binding method [27]. Assuming that 
the important features of the band structure are due to Co and O ion interactions, we used a 
basis set consisting of 5 Co-d orbitals for each of two Co sites and 3 O-p orbitals for each of 
four O sites resulting in 25 Slater Koster (SK) hopping parameters and 5 on-site energies to 
vary. The six bands (4 e g i and 2 ai g ) nearest the Fermi energy are cleanly separated from the 
e g bands above and O 2p bands below by an energy gap [26] allowing them to be isolated from 
all other bands during the fit by a weighting of the eigenvalues. The orthogonality of our basis 
set was enforced and site symmetry taken into account by requiring that on-site energies of 
orbitals in a given irreducible representation be the same. We kept neighbors up to a distance 
greater than half the c-axis lattice parameter to account for inter-planar coupling. Because 
this basis set and number of neighbors results in a large amount of freedom, the parameters 
that fit the LAPW band structure accurately are not unique, and several solutions which 
nicely reproduced the calculated dispersions were found. In order to distinguish the desired 
physically relevant solution (i.e. atomic-like orbital picture), we compared the tight-binding 
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Fig. 1 - A tight-binding fit to the S — 0.7 optimal apical oxygen height band structure. The rectangular 
box indicates the Fermi level crossing responsible for the e g i elliptical hole pockets. The upper inset 
has an expanded c-axis to eliminate inter-planar coupling. 



eigenvectors to the orbitally resolved LAPW band structure at each k-point. The criteria that 
the characters of each eigenvalue must qualitatively match was used to identify a single set of 
parameters as the correct one. 

For each doping, we fit to band structures with three different apical O heights: one above, 
one below and one at the optimal position. Since the NRL-TB parameters are expressed as 
polynomials in distance, we could fit all O positions at each doping simultaneously, thereby 
explicitly including O height dependence in the parameterization. The familiar SK elements 
are then easily obtained for any structure by evaluating the NRL-TB parameters at a specific 
distance. A table of SK elements and on-site energies for all three dopings is provided in 
Tables III and ITT1 with the O position set at its relaxed height in each case. A typical fit (that 
of 2=0.7) is shown in fig. ^ The Fermi level crossings are reproduced almost exactly, though 
there are small discrepancies in dispersion lower in energy, most notably between T and M. 
The rectangular box highlights the e g i band as it breaks above the Fermi energy; this feature 
is responsible for the small hole Fermi sheets which will be the central concern of this paper. 

The TB fit yields relative positions of the on-site energy parameters for the Co-d orbitals 
as expected from simple symmetry considerations: the xy,x 2 -y 2 orbitals are highest in energy, 
the xz,yz orbitals sit lowest and the 3z 2 -r 2 (or ai s ) orbital is raised slightly above the latter. 
Similarly, the in- plane O onsitc energies, p^ and p a , are higher than the out-of-plane energy, 
p z . It is important to note that the orbitals familiarly designated as e s ' and e g are both 
combinations of the xz,yz and xy,x 2 -y 2 representations in the hexagonal coordinate system. 
On a simple triangular lattice (that of the Co ions alone), these representations do not mix, but 
the distorted octahedra of O ions breaks local z-reflection symmetry and facilitates hoppings 
that would otherwise be forbidden, i.e. Co-Co hopping between orbitals that are ostensibly 
orthogonal takes place through O ions positioned non-symmetrically above and below the 
Co plane. A combination of all characters (3z 2 -r 2 , xz,yz, and xy,x 2 -y 2 ) is present anywhere 
that a,\ g and e 9 ' bands mix. At the very top of the band complex, mostly above the Fermi 
energy, the character is purely ai s , but elsewhere throughout the energy region e g > character 
is also present. Both band types cross the Fermi level, necessitating a multi-band treatment 
for accurate reproduction of the DOS and Fermi surfaces. At the very least, any model meant 
to reproduce both a lg and e fl / bands on a triangular lattice must include all five Co d-orbitals, 
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Table I - Hopping Integrals (meV) 



Co-O O-O Co-Co 







dpa 


dpir 


ppa 


Pp-K 


dda 


ddix 


ddS 




1st n.n. 


-1286 


1054 


112 


30 


-326 


34 


34 




2nd n.n. 


-33 


21 


78 


228 


2 


72 


-68 


d=0.3 


3rd n.n. 





-8 


56 


132 





13 


-62 




4th n.n. 






53 


122 










5th n.n. 






18 


30 










1st n.n. 


-1422 


1003 


250 


20 


-405 


63 


36 




2nd n.n. 


-53 


33 


135 


198 


33 


63 


-83 


d=0.5 


3rd n.n. 


-2 


-3 


107 


133 


4 


11 


-57 




4th n.n. 






100 


115 










5th n.n. 






44 


27 










1st n.n. 


-1551 


929 


295 


59 


-433 


90 


16 




2nd n.n. 


-88 


64 


309 


187 


52 


56 


-97 


d=0.7 


3rd n.n. 


-4 


-9 


204 


156 


6 


10 


-47 




4th n.n. 






167 


119 










5th n.n. 






55 


17 









Table II - On-site Energies (eV) 





£=0.3 


:r=0.5 


z=0.7 


■2 '2 

xy,x -y 


3.578 


4.204 


4.762 


xz,yz 


2.422 


3.184 


3.850 


z 2 


2.558 


3.360 


4.122 


Px,Py 


0.272 


0.653 


1.116 




-0.830 


0.0 


0.735 



and the symmetry breaking property of the O ions further implies that even this will not 
necessarily account for all dispersions. 

We now discuss the main doping and related structural dependencies of the band disper- 
sions. In the VCA, the addition of Na narrows the overall bandwidth of what would be, in an 
ideal octahedral environment, the t2g complex. The insensitivity of this bandwidth to c-axis 
compression suggests that extra charge shifts the balance between competing in-plane interac- 
tions, thus decreasing dispersion. Hoppings perpendicular to the plane break the degeneracy 
of the two layers and cause an important splitting of the ai g bands near the T point, but have 
relatively little effect elsewhere. If inter-layer hopping is eliminated by, for example, a pro- 
nounced expansion of the c-axis length, the six bands collapse onto three doubly degenerate 
bands and the splitting disappears (see inset of fig. QJ. Not surprisingly, inter-layer coupling 
increases with growing Na content until at x=0.7, the two ai g bands are separated by as much 
as 20% of the total t2 g bandwidth along the T-A line. The final, and in terms of this discussion 
most important, effect of Na charge is its determination of the apical O height which in turn 
establishes the magnitude of the e g <-ai g split. Coulomb repulsion in the Co-0 plane shifts 
the apical O position upward, easing the octahedral distortion and lowering the ai g band 
weight center relative to e g >. Additional charge mainly enters the lowered ai s band, leaving 
the e g ' hole pockets intact, though diminished in size. The cumulative result of these effects is 
summarized in fig. [3 where the varying bandwidths and inter-planar splittings corresponding 
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Fig. 2 - The bandstructures of Na^CoC^ for three values of x, plotted with Fermi energies aligned. 
The x=0.7 calculation is performed at both relaxed and unrelaxed O positions; in the latter the e g / 
hole surface disappears. 



to each Na level are shown. A second calculation of the x=0.7 compound, but with the apical 
O at an unrelaxed height (we chose the £=0.3 height for comparision) illustrates the effect of 
apical O re-positioning. At the unrelaxed (lower) height, only a lg bands cross the Fermi level 
and e g i states are completely filled, eliminating the small pockets from the Fermi surface. 

We calculated the Fermi surfaces for each doping using our TB parameters and a very 
dense mesh of nearly 8000 fc-points in the BZ. Generating the eigenvalues at each of these 
points in the tight-binding scheme requires less than l/1000th the amount of computing time 
that an LAPW code uses for the same task. For each doping, there is a doubly degenerate 
ai ff -like hexagonal large hole pocket at the top of the BZ, and six doubly degenerate small hole 
pockets near K which are e ff '-like (see fig. |3J). Moving toward the zone center, the two surfaces 
disperse away from each other, one surface growing and the other shrinking. In the :c=0.3 
case, the dispersion is small for both surfaces, and the system is essentially two-dimensional. 
For the £=0.5 doping, the dispersion is stronger, particularly that of the inner or shrinking 
surface which becomes somewhat hourglass shaped. At £=0.7, the z-dispersion of the inner 
surface is so strong that it splits into two separate sheets, disappearing entirely at the zone 
center. The outer surface, in contrast, remains fairly dispersionless. Both surfaces for this 
doping are shown in fig. The increase in z-dispersion with electron number is due to both 
the upward shift of the O ions as more charge is dumped into the Co-0 plane and to increased 
hybridization effects of the charge itself. When Na is added, holes come mainly out of the 
shrinking inner surfaces, mainly the hexagonal one and to a lesser degree the elliptical ones, 
with the effect that the e ff < hole sheets are robust with respect to doping in the outer surface 
and are preserved, though barely, in the inner surface. 

Recent photoemission experiments [28,29] have clearly seen a large Fermi surface consistent 
with the a\g derived section predicted in LDA calculations. The small primarily e' g derived 
pockets were not seen. This could be due to surface sensitivity of these states, for example, 
due to a different surface potential related to surface Na stoichiometry, which would shift the 
positions of the these bands or to a relaxation of the O height in the top layer, which would 
also shift these bands. As we have demonstrated, the apical O must be a certain distance 
from the Co plane for e g < Fermi sheets to exist at the x=0.7 doping. Therefore, contraction 
of the Co-0 complex at the surface may indeed serve to depress the e g > feature responsible 
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Fig. 3 - The outer (left) and inner (right) Fermi sheet of Nao.7Co02. The increased charge present 
in the Co-O planes both raises the apical oxygen position and increases hybridization, making the 
z-dispersion non-negligible. The inner surface shows 3-d dispersion in large (ai 9 ) and small (e g /) hole 
pockets 



for the small hole pockets as exemplified by the unrelaxed calculation of Nao.7Co02 in fig. 121 
Alternatively, it may be that correlation effects neglected by the LDA change the position of 
the e g i band and eliminate the small pockets at dopings where ARPES measurements have 
been made, namely at £=0.65 and x=0.7. K.-W. Lee et al have suggested, based partly on the 
same experimental evidence that prompted division of the phase diagram into three sections, 
that non-negligable static correlation exists only in the high Na region [30]. It is possible 
to produce a relative shift of the a g and e' g derived bands in our tight binding analysis by 
shifting the corresponding on-site parameters for x > 0.5. This mimics what would happen in 
a non-spin-polarized LDA+U calculation [31] where the tendency toward integer occupation of 
orbitals strongly disfavors small Fermi surfaces such as the e s < pockets. We caution that while 
LDA+U band shifts have been used e.g. in Ni to improve the agreement between calculated 
and experimental Fermi surfaces [32], the physics of applying a static mean field correlation 
like LDA+U to a metallic system can be questioned [33] , as it is known that fluctuation effects 
as in the dynamical mean field theory (DMFT) are very important. A 0.067 Ry shift of the ai g 
onsite parameter (corresponding to a Hubbard U of 3 eV) eliminates the small ellipses from the 
2=0.7 Fermi surface and enlarges the hexagonal pocket slightly to accomodate the lost holes. 
If these hole pockets are indeed missing due to correlation effects, then the bandstructure can 
successfully be modelled with a single band (ai g ) and the distinction between triangular and 
pseudo-hexagonal structures discussed previously becomes unimportant. 

In conclusion, we have constructed a tight-binding Hamiltonian which can accurately re- 
produce the LDA Fermi surfaces and band dispersions of Na^CoC^ for x=0.3, 0.5 and 0.7. 
Relaxation of the apical oxygen position at each doping reveals that the Fermi surface is 
non-trivially dependent on doping level and that holes preferentially enter and leave from the 
larger ai 9 -like surface, preserving the smaller hole pockets as charge is added. Foo et al [5] 
observed three rather distinct regions in the phase space of x, two magnetic and one insu- 
lating. Our calculations show that the system passes from quasi-two-dimensionality at low 
Na content to anisotropic three-dimensionality at high Na content. This may have bearing 
on the very different magnetic and thermodynamic properties observed on either side of the 
insulating phase boundary, though the insulating state itself is not reproduced by the LDA. 
We address the differences in predicted and observed Fermi surfaces and postulate that either 
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surface distortion or possibly correlation could account for the discrepancies. 
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